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Fluctuating "order parameter" for a quantum chaotic system 
with partially broken time-reversal symmetry 

S. A. van Langen, P. W. Brouwer, and C. W. J. Beenakker 
Instituut-Lorentz, University of Leiden, P. 0. Box 9506, 2300 RA Leiden, The Netherlands 

The functional p = \ J drip 2 \ 2 of a chaotic wave function ip(r) plays the role of an "order parame- 
ter" for the transition between Hamiltonian ensembles with orthogonal and unitary symmetry. Upon 
breaking time-reversal symmetry, p crosses over from one to zero. We compute the distribution of p 
in the crossover regime and find that it has large fluctuations around the ensemble average. These 
fluctuations imply long-range spatial correlations in ip and non-Gaussian perturbations of eigenval- 
ues, in precise agreement with results by Fal'ko and Efetov and by Taniguchi, Hashimoto, Simons, 
and Altshuler. As a third implication of the order-parameter fluctuations we find correlations in the 
response of an eigenvalue to independent perturbations of the system. 



PACS numbers: 05.45. +b, 24.60. Ky, 42.25.-p, 73.20.Dx. 



Wave functions of billiards with a chaotic classical dy- 
namics have been measured both for classical and 
quantum mechanical waves |^|]. The experiments are 
consistent with a Xp distribution of the squared modulus 
|?/>(r)| 2 of a wave function at point r, the index (3 = 1 or 2 
depending on whether time-reversal symmetry is present 
or completely broken. These two symmetry classes are 
the orthogonal and unitary ensembles of random-matrix 
theory (5) . For a complete description of the experiments 
one also needs to know what spatial correlations exist 
between |i/j(fi)| 2 and IV'C'^)! 2 a t two different points and 
how these correlations are affected by breaking of time- 
reversal symmetry. In the orthogonal and unitary ensem- 
bles it is known that the correlations decay to zero if the 
distance |r 2 — r\| greatly exceeds the wavelength A [||. 

Recently, Fal'ko and Efetov |?J managed to compute 
the two-point distribution ^2(^1,^2) in the crossover 
regime between the orthogonal and unitary ensembles. 
(We abbreviate p, = V|^(ri)| 2 , with V the volume 
of the system.) They found that the two-point distri- 
bution does not factorize into one-point distributions, 
P2(pi,P2) ^ Pi{pi)P\{jp2), even if \r 2 - fi| > A. The ex- 
istence of long-range correlations in a chaotic wave func- 
tion came as a surprise. 

Two years earlier, in an apparently unrelated paper, 
Taniguchi, Hashimoto, Simons, and Altshuler || had 
studied the response of an energy level E{X) to a small 
perturbation of the Hamiltonian (parameterized by the 
variable X). They discovered a non-Gaussian distribu- 
tion of the level "velocity" dE/dX in the orthogonal to 
unitary crossover. This was remarkable, since the distri- 
bution is Gaussian in the orthogonal and unitary ensem- 
bles. 

It is the purpose of the present paper to show that 
these two crossover effects are two different manifesta- 
tions of one fundamental phenomenon, which we identify 
as order parameter fluctuations. The "order parameter" 
is the real number p — \ J dri[i 2 \ 2 in the interval [0,1], 
which equals 1 (0) in the orthogonal (unitary) ensemble. 



The possibility of fluctuations in p was first noticed by 
French, Kota, Pandey, and Tomsovic ||, but the distri- 
bution P{p) was not known. We have computed P{p) in 
the crossover regime, building on work by Sommers and 
Iida and find a broad distribution. Previous the- 
ories for the crossover by Zyczkowski and Lenz EL by 
Kogan and Kaveh and most recently by Kanzieper 
and Freilikher |ll| amount to a neglect of fluctuations in 
p, and thus imply the absence of long-range correlations 
in ip(r) and a Gaussian distribution of dE/dX. Con- 
versely, once the fluctuations of the order parameter are 
properly accounted for, we recover the distant correla- 
tions and non-Gaussian distribution of Refs. |?]||, and 
find a novel correlation between level velocities for inde- 
pendent perturbations of the Hamiltonian. 

We start from the Pandey-Mehta Hamiltonian [p^JT^] 
for a system with partially broken time-reversal symme- 
try, 



H = S + ia(2Ny 1/2 A, 



(1) 



where a is a positive number, and S (A) is a symmetric 
(anti-symmetric) real N x N matrix. The matrix S has 
the Gaussian distribution 



P(S) cx exp (-^Nc^TrSS^ 



(2) 



and the distribution of A is the same. The real parame- 
ter c determines the mean level spacing A at the center 
of the spectrum for A ^S> 1, by c = A A/71- . The distribu- 
tion of H crosses over from the orthogonal to the unitary 
ensemble at a ~ 1. The wave function ipk of the fc-th 



energy level at widely separated points 



fl » A) 



is represented by the unitary matrix U that diagonalizes 
H: 



tk ■ 



(3) 



Consider now an eigenvector \u) = (Uik, U-ik, ■ ■ ■ , U Nk)- 
(Since we deal with a single cigenstate, we suppress the 
level index k.) Following Ref. (^] we decompose \u) in 
the form 



1 



J4> 



(t\R)+i^l-t 2 \I) 



(4) 



where \R) and \I) are real orthonormal TV-component 
vectors, and <f) £ [0, n/2) and t £ [0, 1] are real numbers. 
This decomposition exists for any normalized vector \u) 
and is unique for t 7^ 0,1. The order parameter p is 
related to the parameter t by 



(2t 2 



(5) 



In the orthogonal ensemble t = or 1, hence p = 1, while 
in the unitary ensemble < = \J\j2 hence p = 0. In the 
crossover between these two ensembles the parameter p 
does not take on a single value but fluctuates. 

To compute the distribution P(p) we use a result of 
Sommers and Iida ]lC| ], for the joint probability distribu- 
tion of an eigenvalue E and the corresponding eigenvector 
\u) of the Hamiltonian (jl]). Substitution of the decompo- 
sition ([|), and inclusion of the Jacobian for the change 
of variables from \u) to p, gives the expression 



P(p)<x 



(1 - p) N / 2 ~V 2 
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We have set E = 0, corresponding to the center of the 
spectrum. We still have to take the limit N — > 00. Ex- 
pansion of Z/v(0) in a scries, 



A' 



k=0 



N-k 



1 d k _1 _1 

a fc = MS5f (l- W )- a (l + W )- a 



fc»i 



U1 = 



(7a) 



^ (7b) 

7T 



and replacement of the summation by an integration, 
yields 



„2JV 



Z N (0) 



Y 2 N N-3/2 



2q 



-erf(io) (8) 



(l+w)6 H 



(6a) 

(6b) 
(6c) 
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for AT> 1. Here erf(ia) = 2in- 1 / 2 f£ e y 'dy. The dou- 
ble energy derivative of Z^{E) is computed similarly, but 
turns out to be smaller by a factor N and can thus be ne- 
glected. The derivatives with respect to b± can be found 
by differentiation of Eq. (ph. Collecting all terms, we find 




erf(ia) 



lair? 



erf(ia) 



(9) 



FIG. [j]. Distribution of the order parameter p for 
a = 1/4, 1, and 4, computed from Eq. (^). The crossover 
from the orthogonal to unitary ensemble occurs when 
awl, and is associated with large fluctuations in p 
around its ensemble average. 

It remains to show that the long-range wave-function 
correlations and non-Gaussian level-velocity distribu- 
tions of Refs. follow from the distribution P{p) 
which we have computed. We begin with the wave- 
function correlations, and consider the n-point distribu- 
tion function 



In Fig. |l] the distribution of p is plotted for three val- 
ues of the crossover parameter a. It is very broad for 
a = 1, and narrows to a delta function at 1 (0) for a — > 
(q — > 00). 



Pn(pi,P2,...,Pn) = (j[S(pi-N\U ik \ 2 )\. (10) 

We substitute the decomposition (Q) and do the average 



2 



in two steps: First over \R) and \I), and then over t. Due 
to the invariance of P{H) under orthogonal transforma- 
tions of H, the vectors \R) and \I) can be integrated out 
immediately. In the limit N — > oo, the components of the 
two vectors are 2N independent real Gaussian variables 
with zero mean and variance 1/N. Doing the Gaussian 
integrals we find a generalization of results in Refs. PJTl| 
to n > 1: 

-1 n 

Pn{pi,P2,---,Pn) = dpP(p)T\F(pi,p), (11a) 

J ° i=l 

F(p,p)=(i- P )-hxp(^y (j^Y (lib) 

Here 1$ is a Bessel function. We see that long-range spa- 
tial correlations exist only if the distribution P(p) of p 
has a finite width. For example, the two-point correla- 
tor (pfpf) — (p? )(?>!) equals the variance of p. The ap- 
proximation of Ref. |ll[] (implicit in Refs. |l^,|l3|) was to 
take p fixed at each a. If p is fixed, P n (pi, ■ ■ ■ ,Pn) —> 
P\{pi) ' ' ' PiiPn) factorizes, and hence spatial correla- 
tions are absent. If instead we substitute for P{p) our 
result (ph, we recover exactly the results of Fal'ko and 
Efetov jlll. 

We now turn to the level-velocity distributions. We 
consider perturbations of the Hamiltonian (|l|) by a real 
symmetric (anti-symmetric) matrix S' (A'), 



H' = H + x a S' + x u i 



A'. 



(12) 



Here X u • X o 9X G real infinitesimals, which parameterize, 
respectively, a perturbation that breaks or does not break 
time-reversal symmetry. The corresponding level veloci- 
ties 



dE k 
dx ' 



dEk 



have distributions 



P(v ) = (S(Vo-^2 U *k U jkS'ji 
\ i,3 

p(v u ) = (s(v u -J2u lk u; kt A' ll 



(13) 

(14a) 
(14b) 



We substitute the decomposition (|J) for the eigenvec- 
tor Uik of H and average first over S' and A', assuming 
a Gaussian distribution for these perturbation matrices. 
After averaging over S' and A', the eigenvector enters 
only via the parameter p. One finds 



P(v ) = / dpP{p)G 1+p {v ), 



(15a) 
(15b) 



where G\± p is a Gaussian distribution with zero mean 
and variance 1 ± p. We have normalized the velocities 



P(v u ) = / dpP(p)G 1 . p (v u ), 



such that v% = v% = 1 in the unitary ensemble. Sub- 
stitution of Eq. (H) for P(p) shows that the distribution 
of v Q coincides with the result of Ref. ||. However, our 
P{v u ) is different. This is because we have chosen A and 
A' to be independent random matrices, whereas they are 
identical in Ref. ||. Independent matrices A and A' are 
appropriate for a system with a perturbing magnetic field 
in a random direction; Identical A and A' correspond to 
a system in which only the magnitude but not the di- 
rection of the field is varied. Eq. ([l5]) demonstrates that 
P(v Q ) and P{v u ) are Gaussians in the orthogonal and uni- 
tary ensembles, since then P(p) is a delta function. In 
the crossover regime the distributions are non-Gaussian, 
because of the finite width of P(p). The relation (|TJ) 
between the distributions of v and p for the GOE-GUE 
transition is reminiscent of a relation obtained by Fyo- 
dorov and Mirlin for the metal-insulator transition [|16| . 
The role of the order parameter p is then played by the 
so-called inverse participation ratio / = J dr . In our 
system NI — > p + 2 for N — > oo. A difference with Ref. 
pa is that our perturbation matrices are drawn from 
orthogonally invariant ensembles, whereas their pertur- 
bation is band-diagonal. 

As a final example of the importance of the order- 
parameter fluctuations in the crossover regime, we con- 
sider the response of the system to two or more indepen- 
dent perturbations, 



H' = H 



i=l 3 =1 



■ujiAj. 



(16) 



For example, one may think of the displacement of m 
different scatterers, or the application of a localized mag- 
netic field at n different sites. Proceeding as before, we 
obtain the joint probability distribution of the level ve- 
locities v Q i — dEk/dx i and v U j = dEk/dxuj, 



P(v x,v o2 , ■ ■ ■ ,v om ,v ul ,v u2 , ■ ■ ■ ,v nn ) 

- / d P P(p)l[G 1+p (v ol )l[G^ p (v 

Jo „•_-, 



(17) 



We see that as a result of the finite width of P(p), the 
joint distribution of level velocities does not factorize into 
the individual distributions (plf), implying that the re- 
sponse of an energy level to independent perturbations 
of the Hamiltonian is correlated. 

To summarize, we have computed the distribution of 
an "order parameter" for the breaking of time-reversal 
symmetry in a chaotic system, defined as the squared 
modulus of the spatial average of the wave function 
squared. Fluctuations of the order parameter from one 
wave function to another exist if time-reversal symmetry 
is partially broken. We have shown that these fluctu- 
ations imply long-range wave-function correlations and 
non-Gaussian eigenvalue perturbations, thereby unifying 
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two previously unrelated discoveries . A novel man- 
ifestation of the order-parameter fluctuations is the ex- 
istence of level- velocity correlations for independent per- 
turbations of the system. 

Note added: We have learned that Y. Alhassid, J. N. 
Hormuzdiar, and N. D. Whelan have been working on 
this same problem, with some overlap of results. 

The authors thank Y. Alhassid, K. B. Efetov, V. I. 
Fal'ko, and S. Tomsovic for valuable discussions. This 
research was supported by the Dutch Science Founda- 
tion NWO/FOM. 
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